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文章 编号 :1005-3085(2011)03-0323-12 


带 有 Beddington-DeAngelis 功能 反应 、 脉 冲 、 连 续 时 沸 
和 广义 扩散 函数 的 捕食 者 - 食 乌 系统 的 定性 分 析 * 


= Ri, PR, 0 X 
(陕西 师范 大 学 数学 与 信息 科学 学 院 ， 西 安 710062) 
4 B: 本 文 定性 分 析 了 具有 Beddington-DeAngelis 功能 反应 、 脉 冲 、 连 续 时 滞 和 广义 扩散 函数 的 捕食 
者 - 食 饵 系统 . 利用 脉冲 微分 方程 的 比较 原理 给 出 了 系统 持续 生存 的 条 件 ， 并 使 用 不 动 点 理论 证 
明了 正 周期 解 的 存在 性 ， 进 而 给 出 了 系统 存在 正 周期 解 的 充分 条 件 . 最 后 通过 构造 Lyapunov 泛 
函 证 明了 系统 周期 解 的 全 局 渐 近 稳定 性 ， 该 结论 可 为 现实 的 生物 资源 管理 提供 可 靠 的 策略 依据 . 
关键 词 : 捕食 者 - 食 饵 系统 ， 脉 冲 ， 时 滞 ， 正 周期 解 ， 全 局 渐 近 稳定 
分 类 号 : AMS(2000) 92D30; 34C10; 34K15 中 图 分 类 号 : 0175.7 文献 标识 码 : A 


1 引言 


捕食 者 - 食 饵 系统 是 非常 重要 的 生态 系统 ， 该 系统 正 周期 解 的 存在 性 和 一 致 持久 性 已 有 相 
关 的 研究 成 果 [" 引 ， 在 种 群 的 相互 作用 中 时 灌 是 不 可 避免 的 ， 文 献 [4-13] PIA T EAE 
者 - 食 饵 系统 的 影响 .迁移 现象 是 生物 种 群生 存 过 程 中 一 种 非常 普遍 的 现象 ， 再 考虑 到 人 为 干 
预 种 群生 长 、 繁 衍 的 情况 ， 文 献 [14-20] 中 研究 了 扩散 及 脉冲 因素 对 捕食 者 - 食 饵 系统 的 影响 . 
最 近 的 研究 表明 ， 在 某 些 情况 下 ，Beddington-DeAngelis 功能 反应 函数 可 以 更 好 的 反映 多 个 食 
饵 和 捕食 者 的 捕食 者 - 食 饵 系统 ， 文 献 [21-26] 讨论 了 Beddington-DeAngelis 功能 反应 函数 在 一 
些 生态 模型 中 的 应 用 
本 文 主要 考虑 下 面具 有 Beddington-DeAngelis 功能 反应 、 脉 冲 、 连 续 时 滞 和 广义 扩散 函数 
的 捕食 者 - 食 饵 系统 
tA#tr, k=1,2,---, 
a = 21[a1(t) — aix (t)z (t) — li (t) f°. ku (s)ai(t + s)ds 
~ BOVE BCs) HWE) MM 
Ly, = ti [ai(t) — aü(t)zi(t) — h(t) f?. kyi(s)oi(t + s)ds] 
T di(t)fií(zi,22,:::,24)) t= 2,3,---,n, (1) 
y = y| — bi(t) — bis (t)y(t) — las (t) fA, ks (s)y(t + s)ds + zou esl: 
t=tk, k=1,2,..., 
i(t) = (1 — 0%)ri(te), i=1,2,...,n, 
y(tz) = (1 — pe) y(t), 
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RP z(t), i= 1,2,… ,mg 的 分 别 表示 食 饵 和 捕食 者 的 种 群 密度 ， 而 且 食 饵 zi( 思 和 捕食 
者 y(t) 被 限制 在 斑 块 1 中 ,a;(t), aalt), l(t), i= 1,2,--- ,n, br (t), ba(t), ct), g(t), a(t), (t), 
7( 坟 是 严格 正 的 有 界 连 续 函 数 且 具有 正 周期 w. di(t) 是 扩散 系数 ， 常 数 7 e [0,-+eo)， 函 
Bki(t) > 0, i=1,2,---,n, kalt) > OREN [—7,0] 上 的 分 段 连续 函数 ， 而 且 是 正规 化 的 
函数 ， 即 


0 0 
/ ku(t)dt => 1, i 一 1, 2, dE: ,n, / kz1(t)dt 一 1. 
广义 扩散 函数 f(z1, zz，… , tn) (i = 1,2, --- ,n) YE [0, +00] 上 连续 且 满 足 如 下 条 件 : 


1) EF > Tj (j Ft, j = L2 ,n) 时 ， fií(z1,22,:-- ; Tn) « 0, fj(z1,22,:*- , Tn) > 0; 
2) H ti < Tj (j #14, j =1,2,--- n) 时 ， fi(z1,22,7-- à) > 0, f;(zx1,22,77- ,Tn) <0; 


3) |fi(zu 22, En) — fií(£5 82, ,£)| < Li Yo lc; zi HPL, ABR. 
Si 


wz 

SRO, i — 1,2, ,n, EES, Hefe ERE, MEO, = 0 i= 

1,2,:,n, Hkkq = Kkytktg = tk tw, KEN. 根据 实际 意义 可 知 , 1—01 > 0, i = 
12,- jn, 1— uk > 0. 参数 如 满足 


O<ti<tea<---<t,<-++-, lim tk =+00. 
大 一 十 co 


H (@- 6%), i=1,2,...,n, , Tl (1 — 44) 是 以 w 为 周期 的 周期 函数 ， 对 所 有 t > 0， 存 
<tk <t 


Oct, «ct 


在 两 个 正常 数 nil 和 Ni WE 


nı S II (1-64) <N, i=1,2, n. 
0<tk<t 


对 所 有 t > 0, 存在 两 个 正常 数 ns。 和 No, WE no < I] (1 — px) € No. id 


Octy «t 
BR 一 [o 22, -- Zn. y) € R^" | > 0, i= 1,2,--: 0. y 2 O}. 
系统 (1) 的 初始 条 件 由 以 下 函数 给 出 
Ti(s) = pi(s) >0, i=1,2,..,n, y(s-—v(s) >0, se[-r,0], (2) 
其 中 Qi, Y € C!([—r, 0], R4+ 1 )， i= 1,2,- ,nN 
为 了 讨论 方便 ， 采 用 以 下 记号 和 假设 . 记 PC(R+,R) 是 满足 以 下 条 件 的 函数 集合 p : 
R* 一 R. Kot) itte R* Htz t, MEM, At, e Rt 是 函数 的 第 一 类 不 连续 点 且 在 该 
点 处 的 左 极限 存在 . 记 带 有 正 周 期 w AY Banach 空间 为 
PC, = (6 € PC((0,u],R) | $(0) = $(w)} x [lléllec = sup 1901}, 
tE[0,w] 
id 
f= min fe), f= max 0, Pel[ goa. 
te[0,u] t€[0,c] w Jo 
其 中 f(t) e PC,. 
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2 系统 (1) 的 一 致 持久 性 


定义 2.1 如 果 存 在 有 界 紧 集 Te R4+1+， 使 得 系统 (1) 满足 初始 条 件 (2) 的 每 一 个 解 最 终 都 
进入 并 滞留 在 集合 中 ， 则 称 系统 (1) 是 一 臻 持久 的 . 
引 理 2.1 假设 有 如 下 的 脉冲 方程 


s'(t) = s(t)(as(t) — ax1(t)s(t)), ttk, k=1,2,---, T 
s(t?) =mas(t), t=th, k-12,-, 
和 非 脉冲 方程 
z' (t) s x(t) (att) — a11(t) Il mya(t)), (4) 


Oxty «t 
其 中 函数 a1(t), a(t) 是 严格 正 的 有 界 连续 函数 且 具 有 正 周 期 w， 参 数 妇 满足 


0< 有 < 如 < .< 大 人:……， lim tx = +00. 
天 一 十 co 


脉冲 参数 mx 满足 0 < mk < 1， 且 存在 一 个 正 整 数 q， A mia = Mk, tptg = tk +w, k € N 
立 ， 则 有 以 下 结论 ; 
1) WR a(t) Æ (4) Æ [—7, +00) EMF, W s(t) = dl RENAE (3) 的 解 


2) WR s(t) Æ (3) Æ [-r, +00) EKAR, Meos TI LEA 的 解 . 
€ty «t 
证 明 设 z(t) 是 (4) 在 [-7, 十 co) 上 的 解 ， 则 s(t) =  [[ mex(t) ERM (tz, tk -- 1], k = 
O<t,<t 
… PIES, HWER ttr k=1,2,---, A 


s'(t) — s(t) (a1 — a11s(t)) = ll Mk r — x(t) (a 一 all I] mxz(t)) ) =0, (5) 


0<tk<t 0<tk<t 
另 一 方面 ， 对 任意 如 € (t, k=1,2 p fü 
s(tg) = Jim. II m= [| melts), 


tk O<t.<t O<t;<tk 


由 已 知 的 脉冲 点 的 条 件 (tt) = I mjz(t). RAD ERATE 


O<t; <tr 
s(tz) = m«s(ta), (6) 

从 (5) 和 (6) 可 知 ， o(t) = JT Tux (t) E (3) 的 解 ， 
x oe dar ec tx em. WU s(t) ZEDX IRI (ty, tei] (5 = 1,2,… ) 是 连续 的 .结合 


z(tt)= JI stt) = II — s(t) = (te), 


Oct; Stk J 0xt; «t. J 


s(t)- [I Zst) = Il s(tk) = (te), 


O<tj<t. 7 Oct;ct, 7 
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MAT a(t) Æ [—7, +00) 上 连续 ， 且 满足 方程 (4)， 所 以 z(t) 


引 理 2.2 Rit) 是 系统 (1) 的 不 变 集 . 

引 理 2.3 设 常 数 M > 0, MABE (1) 满足 初始 条 件 (2) 的 任 一 解 (ms (£), x2(3),… ,zn(t), y(t))， 
当时 间 t 充 分 大 时 有 zi(t) < M, i=1,2,:…,n, y(t) € M. 

证 明 令 V(t) = max{zx1(t),7z2(t),… ,zn(t)}， 沿 系统 (1) 计 算 V(t) 的 右上 导数 ， 则 有 以 下 
几 种 可 能 : 

1) V(t) 2zi(t), Bu 


ID 过 s(t) 是 (4) 的 解 . 


0<t 


D*V(t) < ax(t)(ax(t) — oan(t)z(t)) € V (E) (a? — a, V (t); (7) 
2) V(t) = z;(t, j—2,3, -,n. W 
DHV (t) € z;(t)(a;(t) — a5; (t)z(t)) < VE) (or — a5; V(t)). (8) 
& Ay € (a, i=1,2 ,n), Án € (d i 1,2, ,n}， 而 且 满 足 
ES 
Ai 
则 结合 系统 (1) I (7), (8) 式 ， 可 以 得 到 如 下 脉冲 系统 


am 
= max { H, i=1,2,---,nh=M, 
ay 


| D*V(t) <V(t)(A.-AnV@), t#te, k-L2.-, m 
V(t*) = mV (t), t=te, k=1,2,..., 
其 中 mx = max{1 — 0i, i =1,2, ,n}， 初 值 条 件 为 
V(s) = max{pi(s), i=1,2,---,n, s €[-7,0]}, , V(0*) = max{y,(0), i = 1,2,.… ,n}. 
考虑 如 下 非 脉冲 方程 
v'(t) = v(t) (A -ån [I mxv(t) ), (10) 


0<tk<t 
其 初始 条 件 为 v(0) = V(0). 对 (10) 式 再 利用 脉冲 参数 的 已 知 条 件 进 行 放大 ， 可 得 
v'(t) < v(t) (Ai ET Ayinv(t)). 


经 过 计算 可 知 ， 存 在 五 > 0, St > TN, v(t) € Mi/m,. BAK HR eS 及 引 
38 2.1 可知 ， “4t>Ty RY, V(t) x NAMi/n; BI25 t Tj 时 ， 有 


zi(t) < NMi/ni1, 4=1,2,--- ,n. (11) 
从 系统 (1) 的 第 三 个 等 式 及 已 知 条 件 ， 可 以 推出 
yO € yO — 50 — 6090 + AS) «vo (S - 4 - tato). (12) 


g” — bi! 
B'bj, 





Bi = ri — b}, By =b), Mə = 
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则 结合 系统 (1) 和 (12) 式 ， 可 以 得 到 如 下 脉冲 系统 


ie cp n k212,-, T 
y(t*) =(1—pe)y(t), t=te, k=1,2,---, 
初 值 条 件 为 y(s) = v(s) > 0, s € [-7,0], y(0*) = v(0). 
考虑 如 下 非 脉 冲 方程 
v'(t) = v(t) (5i -Bi J] a- inyv(t)), (14) 


O<t,<t 
其 初始 条 件 为 v(0) = v(0). 对 (14) 式 再 利用 脉冲 参数 的 已 知 条 件 进 行 放大 ， 可 得 
v'(t) < v(t) (Bi = Byngv(t)). 


xiu A. FET > Ti, 4t > Toit, v(t) € Meo/no. 由 脉冲 比较 定理 H6] 及 引 
理 2.1 可知 ， 当 t+ > 时， 有 





y(t) < N2 M2 /n2, (15) 
4 M = max{NiMif/ni,NoM2/n2}， 则 从 (11) 和 (15) 式 可 得 : 当 上 > T2 B, A 
zit) <M, i=1,2,:...,n, wt)zM, (16) 
即 当 t 充 分 大 时 ,zi(t) € M; i=1,2,---,n, y(t) x M. 
l I al 
— MIT 一 cm . — MI® 
m -(2 Mii 5 Hr d m li i=2,3, 5n] 
a Qj 
定理 2.1 #8 
l 
(H1) mı > 0; (H2) DN op P nin TNM — bT — MIF; > 0; 


则 系统 (1) 是 一 致 持久 的 . 
证 明 W (x1 (t), z(t), ms (0, (0) 是 系统 (1) 满足 初始 条 件 (2) 的 解 ， 则 从 系统 (1) 的 第 
一 个 和 第 二 个 等 式 和 引 理 2.3， 可 知 存在 Ts > 0, Ut» TEL. A 


z, > xila} — MI — c" 4l — az) 4- di fi(21, 22, Zn); 
a > xi (al pem Mif 一 a; xi(t)) + dlfi(zi 22, vA Za), i= 2,3, UO, 
WV, (t) = min(zi(t), z2(t),--- ,zn(t)}， 沿 系统 (1) 计 算 Vi(t) 的 左下 导数 ， 则 有 以 下 几 种 
可 能 : 
1) Vi(t)=ai(t), WW 
D-Va(t) > Vi (t) (ay — MIT - e" /5! — aV (t); (17) 
2 V(t) = z(t), i=2,3,---,n, 则 
D-Vi(t) 2 Vi (ai - MUP — ag Vi(t)). (8) 
4 à; € {a} - MIM — c" /4l, al — MIB, i = 2,3,--- n}, õn € (a, 1=1,2,---,n}, 而 且 
满足 
a , (s - MIR —c™/ adl - MI i 


1 
=— = min , 


O11 anu ay 


=1,2,--- in} =m, 
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则 结合 系统 (1) 和 (17), (18) 式 ， 可 以 得 到 如 下 脉冲 系统 
D-Vi(t)>V()(Gi-GuVit)), tAte, k=1,2,---, 
| Vi(tt)=pVi(t), t=te, k=1,2,---, 
其 中 px = min{1 — 64,1 =1,2,--- ,n}. 初 值 条 件 为 
Vi(s) = min(ei(s), i = 1,2,--- se [-7,0]}, V(0*) = min{y;(0), i =1,2, , rn}. 
考虑 如 下 非 脉冲 方程 


(19) 


vi (t) = vi(t) (a — G11 II pevi(t)), (20) 
O<t,<t 
其 初始 条 件 为 v1(0) = Vi(0). 对 (19) 式 再 利用 脉冲 参数 的 已 知 条 件 进行 缩小 ， 可 得 
vi(t) 2 vi (£) (à — dNiv1(t)). (21) 


经 过 计算 可 知 ， 存 在 到 > TS, Mit Tah, v(t) > m1/(2N1)， 由 脉冲 比较 定理 (19 及 引 
理 2.1 可 知 ， EE; > T4 时 ， Vi(t) > nym,/(2N1)> BI 34 t > Ta 时 ， 有 


nym 
(> 
zi(t) = 2N, , 


从 系统 (1) 的 第 三 个 等 式 ， 引 理 2.3 和 (22) 式 ， 可 知 当 t > UN, A 


l 
(ag umo 
y (t) > y(t) us + B™nym, + 2MNiy™ 





i21,2,--,n. (22) 


- by - Mig - wt). (23) 
* 


= l E b 
SE — bY — MI, bj = bii» m2 = - , 
11 


rs 2Nia™ + B™nym, + 297M M 
则 结合 系统 (1) 和 (23) 式 ， 可 以 得 到 如 下 脉冲 系统 


| vt) > wt)(bh-buwt)) tte k-212;-, i5 
yt*) = (1— pr)y(t), t=th, k=1,2,---, 
其 初 值 条 件 为 y(s) = v(s), s € [-7,0], y(0*) = v(0). 
考虑 如 下 非 脉 冲 方程 
w(t) = ult) (b -bu I[ G- u), (25) 


O<tk <t 


其 初始 条 件 为 w(0) = v(0). 对 (25) 式 再 利用 脉冲 参数 的 已 知 条 件 进 行 缩小 ， 可 得 
u'(t) > u(t) (fr — b11 Nou(t)). (26) 


经 过 计算 可 知 ， 存 在 站 > T4. St > Ts, ult) > ma/(2N3). 由 脉冲 比较 定理 06} 及 引 
理 2.1 可 知 ， 当 上 > TS B, A 


nem 


> x 
y(t) > 2N; 





(27) 
4 m = min{ H oN, k ON: nema}, WI 
T = {(a1(t), zo(t),--- ,an(t), y(t)) |m < zit) € M, i=1,2,--- ,n, m < y(t )<M}. (28) 


BARATERLT 的 一 个 有 界 紧 子 集 . 从 以 上 讨论 可 知 ， 当 + 充分 大 时 系统 (1) 满足 初始 
条 件 (2) 的 每 一 个 解 最 终 都 进入 并 滞留 在 集合 工 中 ， 所 以 系统 (1) 是 一 致 持久 的 . 
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3 ” 正 周 期 解 的 存在 性 及 稳定 性 


本 节 中 首先 给 出 系统 (1) 存在 正 周 期 解 的 充分 条 件 ， 然 后 讨论 周期 解 的 稳定 性 . 

定理 3.1 如 果 条 件 (H1) 和 (H2) 成立， 那么 以 w > 0 为 周期 的 系统 (1) 至 少 存在 一 个 正 周 
期 解 . 

证 明 由 引 理 2.2 可 知 ，R?t! 是 系统 (1) 的 不 变 集 ， 因 而 系统 存在 一 个 正 不 变 集 ， 则 
A (28) 中 所 定义 的 集合 工 是 系统 的 一 个 正 不 变 集 . 

定义 Poinare 映 秧 B:T — D, Bp 


$ (zi (0*), zo (0*), db ,Zn(01), y(0*)) m (zı(wt), za(wt), wet :In(wt),y(w*)), (29) 


MAREA DIÉRAIU LHARAGCHATR, SUN ORTSIT MN ABU. 由 解 对 初 值 的 连续 依赖 
性 ， 下 为 连续 算 子 ， 由 Brouwer 不 动 点 定理 可 知 ， 更 在 工 中 存在 不 动 点 ， 即 系统 (1) 至 少 存在 
一 个 正 周 期 解 . 

定理 3.2 (F(t), Z2(t),--- , 24 (0), H(t) 是 系统 (1) 满足 初 值 条 件 (2) 的 一 个 正 的 有 界 解 . 
车 有 条 件 (H1), (H2) 及 以 下 条 件 成 立 


c™B™M + go + g^" M n qn; 
j-2 
D. d"L; 
(H4) ah- iR- 9$, — >0, i=2,3,---,n: 
j-ljzi 
(H5) bii +% -— A >0; 


EP Am) = (a! + g'm + ym), WI (Z(t), £2(t),--- 34 (0), JA) 是 全 局 渐 近 稳定 的 . 

证 明 WE (a(t), r(t) , an(t), y(t) BABE (1) 满足 初 值 条 件 (2) 的 一 个 解 . 从 定理 2.1 和 
(28) 式 中 定义 的 集合 工 可 知 ， 存 在 T6 > 0, 4t>Te WN, (21 (t), z2(0,--- ,Zn(t),y(t)) € T. F 
在 to > 0， 当 上 > Tg -- to b. (Z(t), 2(t),… 4 (0), 9(t)) Er. 即 当 上 > T6 十 to 时 ， 有 


令 
z; = lnz;(t), 1=1,2,---,n, y* =Iny(t), žy =lnē;(t), i-1,2,..,n, Y= 1ng(t). 


考虑 如 下 的 Lyapunov 函数 
n 0 t 
vO = Y (leto - | +n 小 kid») 人 lxs(v) ~ &(»)dvds) 


i=1 


+ 





0 t 
v0 -g 0m 人 kaa (5) J WO) ~ Rlavas, (30) 
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JJ tE (th-1, th] C [Te + to, to) iF, HARA nA 的 右上 导数 ， 可 得 


0 
D*Va(t) = sign(zi — Zt) | — axi(£)(z1 — £1) — l1 (t) J kii (s) (zi (t + 5) — žı (t + s))ds 





ES c(t)y(t) c(t)g(t) 
alt) + Bait) (D ^ alt) + POE) + IO) 
+O porta En) 一 a PS 4) 





十 S aui — di) | —aü(t)(z; — ZF) — l(t) [i kii(s) (zi(t + s) — ji(t + s))ds 


i-2 


at ) di I 





VT (ri za Bn) — E fias En | 


ssign(y" - j| - hut 9 - tato f. ntt 9) e+ s))as 


LLL gni) : sal) a] 
a) + Be) FOO — a(0) FBO) OO 


n 0 0 
+ m/v) -sls f futt s) = zil pee] 
i=l 2 FN 


0 0 
Haf olio - stolas — f kallu s) - ate as]. (31) 
注意 到 
oao dO OHO 
sien(aj - a) — De pac 359503559 


有 c" am +c? gm M 


mae y(t) oe ee 


ne paccm M ale )- &i(t)], (32) 


: * 一 水 g(t)zi(t) g(t)žı (t) 
sign(y Psy ono rut UN OES OENORSTOTO! 


EE OS 0 + be - a0), (33) 








其 中 A(m) = (a! + Bim + ym). 
4 
-— 








fi(ii, 2, -- i), 


d;(t) 
(£1, T2,' (448) E [3 


Di(t) = sign(x} — à) [ŻE 
从 广义 扩散 函数 fi(z1,22, 7 tn) 的 定义 出 发 ， 分 情况 讨论 xz; 和 ;的 大 小 关系 ， 可 推导 出 


D,(t) < E TEY e - &;. (34) 
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参照 (32)-(34)， 对 (31) 式 进 行 放 大 ， 可 得 


m m gm M 十 gmam + orymM n dm 
Dtva(t) < -(a m TATEM 12s - | 





j=2 





-M(u-m- Y; Te-a 


i-2 jfi 
l 1 ¢™B™M + g"o" + g™y"™M E 
- (bis +l 一 Am c9 ed (35) 


根据 条 件 (H3)-(H5)， 可 知 存在 一 个 常数 5 > 0, 2t 6 (te_1, te] C [Te + to, too) FT, A 


D*V4(t) < -ô X le: - žil + ly - al), (36) 


n 0 £t 
Ve) = Y» (let) - ated nf CON [z(») — i (»)|dvds) 


w= 


0 tt 
been - sep Hm f kao f lue) - stolevas 
-T tk +s 


E ( [inc i 0i) Ti (tk) 一 In(1 uro ki (tx)] 


i=l 


-m ET 人 
tls 


tk 
+% T Kits) f ly(v) — §(v)|duds = lim V(t) = Va(ta). (37) 
tis tot 
由 (36) 式 和 (37) 式 ， 可 知 当 t+ € [Tg 十 如 ,+co) 时 ， 有 
D*Va(t) < (Y le: - a1 ly — il). (38) 
i=1 


对 (38) REKE) [Tg + to, t] LR, WE 


Wes ( 3 [zi(s) — ži(8)| + |y(s) — G(s)| )ds < V2(Te + to) < +00, 
于 是 . 
上 (> leis) — «(9| + ule) — Hs)| )às < V (Ta + to)/6 < +00, 
6+t ` i—1 
所 以 


lzi(t) — &;(t) 十 jut) = g(t)| € Li(Ts + to, +00). 
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由 定理 2.1 可知， 函数 |zi(t) — z;(t), i= 1,2,--- ,n, y(t) — g(£)| YE [To + to, too) LW SK 
有 界 ， 因 此 


D le - i| + yO - 90) 
££ [Te + to, +00) 上 是 一 致 连续 的 ， 
由 文献 25] 可 知 


Y ast) — 2| + ly) - GO| = 0. 
i=1 
由 以 上 讨论 可 知 


im jelt) - z4()| 20, 1=1,2,---,n, ,lm ly) - 9(t)] = 0, 
由 文献 [22] 中 全 局 渐 近 稳定 性 的 定义 ， 可 知 (21(£), 22(0),--- , 24 (0), W(t) 是 全 局 渐 近 稳定 的 . 

使 用 定理 3.1 和 构造 与 定理 3.2 中 相似 的 Lyapunov 函数 ， 可 以 得 到 如 下 推论 . 

推论 3.1 如 果 条 件 (H1)-(H5) 成 立 ， 则 以 w > 0 为 周期 的 系统 (1) 存在 一 个 周期 为 w > 0 正 
周期 解 ， 而 且 是 全 局 渐 近 稳定 的 . 

参照 文献 [18-21] 和 文献 [24] 中 对 捕食 者 - 食 饵 系统 存在 正 周 期 解 的 充分 条 件 ， 不 难 获得 以 下 
关于 系统 (1) 至 少 存在 一 个 周期 为 w > 0 的 正 周 期 解 的 充分 条 件 . 

定理 3.3 ”如 果 以 下 条 件 


4 wor iam 一 
(H6) a 135 m(1 + 6h) > c/y+ Lid: 
k=1 
q : 一 
(H7) t+ $5 n(14 64) > Lidi, i= 2,3,---,n: 
k=1 


(H8) 1 55 In(1 + px) > bis 
k=1 
成 立 ， 则 系统 (1) 至 少 有 一 个 周期 为 w > 0 的 正 周期 解 . 


4 结论 


本 文 研究 了 具有 Beddington-DeAngelis 功能 反应 、 脉 溃 、 连 续 时 滞 和 广义 扩散 函数 的 捕食 
者 - 食 饵 系统 的 一 致 持久 性 和 周期 解 ， 首先 引入 三 个 引 理 ， 然 后 利用 脉冲 微分 方程 的 比较 原理 
讨论 了 系统 (1) 持续 生存 的 条 件 ， 使 用 Brower 不 动 点 理论 证 明了 正 周 期 解 的 存在 性 ， 进 而 给 出 
系统 (1) 存在 以 w > 0 为 周期 的 正 周期 解 的 充分 条 件 . 通过 构造 Lyapnov 函数 证 明了 系统 (1) 的 
周期 解 是 全 局 渐 近 稳定 的 .从 研究 中 发 现 ， 在 满足 一 定 的 条 件 下 脉冲 参数 影响 了 系统 各 种 群 数 
量 的 上 下 界 及 全 局 渐 近 稳定 性 ， 扩 散 行为 不 影响 系统 的 一 致 持久 性 ， 但 在 系统 的 全 局 渐 近 稳定 
性 中 发 挥 了 作用 ， 时 滞 参 数 对 系统 的 影响 不 显著 . 


参考 文献 : 


[1] Cushing J M. Periodic time-dependent predator-prey system[J]. SIAM Journal on Applied Mathematics, 
1977, 32(1): 82-95 
[2] Levin S A. Dispersion and population interaction[J]. The American Naturalist, 1974, 108(960): 207-228 


第 3 期 王 烈 ， 等 : 带 有 Beddington-DeAngelis 功能 反应 和 广义 扩散 函数 的 捕食 者 - 食 饵 系统 的 定性 分 析 333 


[3] 
[4] 


[5] 


[14] 


[15] 


[16] 


[17] 


18] 


[19] 





[20] 
[21] 
[22] 


[23] 


Gaines R E, Mawhin J L. Coincidence Degree and Nonlinear Differential Equations[M]. Berlin: Springer- 
Verlag, 1977 

Cai L M, Song X Y, Chen Q J. Permanence and stability in a delayed ration-dependent predator-prey 
system with stage structure for predator[J]. Chinese Journal of Engineering Mathematics, 2006, 23(3): 
537-542 

SM, DRLE. 四 种 群 食物 链 方程 的 整体 渐 近 稳定 性 [ 吕 . 工程 数学 学 报 ，2006, 23(3): 407-413 

Li Y L, Ma Y C. Global asymptotic stability of four-species food-chain systems{J]. Chinese Journal of 
Engineering Mathematics, 2006, 23(3): 407-413 

秦 发 金 . 一 类 具有 收获 率 和 比率 的 时 滞 阶 段 结构 的 扩散 捕食 系统 的 多 重 正 周期 解 [J]. 工程 数学 学 报 ，2009, 26(4): 
671-679 

Qin F J. Multiple periodic solution for a delayed stage-structure predator-prey systems with ratio- 
dependence and harvestion rate and diffusion[J]. Chinese Journal of Engineering Mathematics, 2009, 26(4): 
671-679 

Yang K. Delay Differential Equations with Applications in Population Dynamics[M]. New York: Academic 
Press, 1993 

Panetta J C. A mathematical model of periodically pulsed chemotherapy: tumor recurrence and metastasis 
in à competition environment[J]. Bulletin of Mathematical Biology, 1996, 58(3): 425-447 

Song X, Chen L. Persistence and global stability for nonautonomous predator-prey system with diffusion 
and time-delay[J]. Computers & Mathematics with Applications, 1998, 35(6): 33-40 

Cuevas C, Pinto M. Asymptotic properties of solutions to nonautonomous Volterra difference systems with 
infinite delay[J]. Computers & Mathematics with Applications, 2001, 42(3-5): 671-685 

Wang L L, Wan T. Existence’and global stability of positive periodic solutions of a predator-prey system 
with delays[J]. Applied Mathématiés and Computation, 2003, 146(1): 167-185 

Ding X H, Lu C. Existence of positive periodic solution of ratio-dependent delay N-species difference 
system[J]. Applied Mathematical Modelling, 2009, 33(6): 2748-2756 

王 爱 丽 ， 陈 斯 养 ， 王 东 保 . D^ X Logistic HARM HEARD THEE). 兰州 大 学 学 报 (自然 科学 版 )，2004， 
40(2): 8-12 

Wang A L, Chen S Y, Wang D B. Asymptotic properties for the general delay Logistic single species 
biological models[J]. Journal of Lanzhou University (Natural Science Edition), 2004, 40(2): 8-12 

王 爱 丽 ， 陈 斯 养 . 具有 多 时 滞 和 广义 扩散 的 捕食 链 模型 的 一 致 持久 性 [可 ]. 南昌 大 学 学 报 (理科 版 )，2006,，30(1): 
17-22 

Wang A L, Chen S Y. Uniform persistence for nonautonomous food-chain system with delays and general 
diffusion[J]. Journal of Nanchang University (Natural Science), 2006, 30(1): 17-22 

李 小 玲 ， 胡 广 平 . IF R-i a RATA EXIT] eA A PS [J]. SINK RER). 2009, 45(2): 
81-84 

Li X L, Hu G P. Stability and periodic solutions in a predator-prey system with delays[J]. Journal of 
Lanzhou University (Natural Science Edition), 2009, 45(2): 81-84 

Bainov D D, Simeonov P S. Impulsive Differential Equations: Periodic Solutions and Applications[M]. New 
York: Longman, 1993 

Xu R, Chen L S. Persistence and global stability for a delayed nonautonomous predator-prey system without 
dominating instantaneous negative feedback[J]. Journal of Mathematical Analysis and Applications, 2001, 
262(1): 50-61 

Tang S, Chen L. The periodic predator-prey Lotka-Volterra model with impulsive effect[J]. Journal Me- 
chanics in Medicine and Biology, 2002, 2(3): 267-296 

Liu X, Chen L. Complex dynamics of holling type II Lotka-Volterra predator-prey system with impulsive 
perturbations on predator[J]. Chaos, Solitons & Fractals, 2003, 16: 311-320 

Cai L M, Li X Z. Positive periodic solutions for à nonautonomous delayed predator-prey system with 
diffusion and impulses[J]. Acta Analysis Functionlis Applicata, 2008, 10(2): 139-149 

Zhang S, Chen L. A study of predator-prey models with the Beddington-DeAnglis functional response and 
impulsive effect[J]. Chaos, Solitons & Fractals, 2006, 27(1): 237-248 

Naji R K, Balasim A T. Dynamical behavior of a three species food chain model with Eddington-DeAngelis 
functional response[J]. Chaos, Solitons & Fractals, 2007, 32(5): 1853-1866 

Huo H F, Li W T, Nieto J J. Periodic solutions of delayed predator-prey model with the Benddington- 


334 T EF Kk 学 学 S 第 28 卷 








DeAngelis functional response[J]. Chaos, Solitons & Fractals, 2007, 33(2): 505-512 

[24] Cai L M, Li X Z, Yu J Y, et al. Dynamics of a nonautonomous predator-prey dispersion-delay system with 
Beddington-DeAngelis function response[J]. Chaos, Solitons & Fractals, 2009, 40(4): 2064-2075 

[25] Beretta E, Takeuchi Y. Global asymptotic stability of Lotka- Volterra diffusion models with continuous time 
delays[J]. SIAM Journal on Applied Mathematics, 1988, 48(3): 627-651 

[26] 王 烈 ， 陈 斯 养 ， 石 茂 . 带 有 脉冲 、 时 滞 和 广义 扩散 函数 的 捕食 者 - 食 饵 系统 正 周 期 解 的 存在 性 [可 . 兰州 大 学 学 报 ( 自 
然 科 学 版 )，2010, 46(1): 96-104 
Wang L, Chen $ Y, Shi M. Existence of positive periodic solution of predator-prey system with impulsive, 
time delay and general diffusion[J]. Journal of Lanzhou University (Natural Science Edition), 2010, 46(1): 
96-104 





Qualitative Analysis of Predator-prey System with 
Beddington-DeAngelis Functional Response, Impulsive, Continuous 
Delay and General Diffusion 


WANG Lie, CHEN Si-yang, SHI Mao 


(College of Mathematics and Information Science, Shaanxi Normal University, Xi'an 710062) 


Abstract: A nonautonomous predator-prey model consisting of n-competing preys and one predator 
with the Beddington-DeAngelis functional response, impulsive, continuous delay and general diffusion 
is proposed. First, it is proved that the system is uniform persistence by using the comparing theorem 
of the impulsive system. Secondly, the existence of periodic solutions is proved through the Brower fixed 
point theory. Through constructing a Lyapunov mapping, the sufficient conditions for the existence of 
the positive periodic solution and the global asymptotic stability of the positive periodic solution are 
obtained. Our results provide a reliable tactic basis for the practical biological resource management. 
Keywords: predator-prey system; impulsive; time delay; positive periodic solution; global asymptotic 
stability 
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